Geometric optics of Bloch waves in a chiral and dissipative medium 
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We present a geometric optics theory for the transport of quantum particles (or classical waves) 
in a chiral and dissipative periodic crystal subject to slowly varying perturbations in space and 
time. Taking account of some properties of particles and media neglected in previous theory, we 
find important additional terms in the equations of motion of particles. The (energy) current density 
field, which traces the geometric optics rays, is not only governed by the Bloch band energy dispersion 
but also involves there additional fields. These are the angular momentum of the particle, the 
dissipation dipole density, and various geometric gauge fields in the extended phase space spanned 
by space-time and its reciprocal, momentum and frequency. For simplicity, the theory is presented 
using light propagation in photonic crystals. 



PACS numbers: 03.65.Sq, 42.25Bs, 03.65.Vf, 72.10.-d 



I. INTRODUCTION 

The transport property of quantum particles (or clas- 
sical waves) in media is one of the central problems in 
physics. A textbook example is the dynamics of Bloch 
electrons in solid state crystals For a long time, the 

semiclassical equations of motion (EOM) of electrons are 
known to be solely determined by the energy band dis- 
persion Recently it was found that geometric gauge 
fields due to Berry phase [3] effects may modify the EOM 
of particles d i] in a chiral medium, leading to new and 
important physical phenomena and applications 0413 ■ 
Here "chiral" means that the time reversal or spatial in- 
version symmetry in the medium is broken. 

However, the EOM in current literature neglect several 
important properties of particles (or waves) and media 
{e.g., self rotation of particles, frequency (or energy) de- 
pendence and dissipation of the media). In this paper, 
we derive a new set of EOM of particles where the ef- 
fects of these properties are included. For simplicity of 
the presentation, we consider a concrete physical system: 
light propagation in photonic crystals although the 
results are applicable to many other systems [ll| . Instead 
of wavepacket dynamics commonly used in previous lit- 
erature, we adopt a geometric optics approach that has 
been proved to be extremely valuable for describing light 
propagation in usual optical materials. However, the tra- 
ditional approach of geometric optics is inadequate be- 
cause the light wavelength is typically the same order as 
the unit cell size of the photonic crystals within which 
the dielectric constant varies. 

The new geometric optics theory is presented for light 
propagation in a frequency band of a chiral and dissi- 
pative photonic crystal under perturbations which vary 
slowly in space and time compared to the lattice con- 
stant and the band gaps. We focus on the Poynting 
vector field, averaged over the crystal unit cell, which 
traces the geometric optics rays {i.e., the EOM). It in- 
volves the photon wavepacket center velocity, and also 
contains two additional terms due to the photon orbital 
angular momentum (POAM) and the dissipation of the 
medium that are absent in previous wavepacket theory of 



photons We show that, in a concrete example, the 

term from POAM can dramatically modify the transport 
properties of light in photonic crystals. The propagation 
of light is also affected by various geometric gauge fields 
due to Berry phase effects. Both POAM and Berry phase 
originate from the chirality of the medium or the parti- 
cle. Two such gauge fields mimic closely the electric and 
magnetic fields for electrons and two other fields 

arise from Berry curvatures in the momentum-frequency 
space, previously called reciprocal electromagnetic fields 
in the context of electron semiclassical dynamics [1, [3] . 
The physical origin of the reciprocal electric field is il- 
lustrated. There are also fields from Berry curvatures in 
other facets of the extended phase space. All these fields 
can be calculated based on the band structure wave func- 
tions of the photonic crystal. 

The paper is organized as follows: Section II intro- 
duces the wave equations and the wavepacket dynamics 
for light propagation in photonic crystals. In section HI, 
a geometric optics theory for light propagation in pho- 
tonic crystals is developed. Section IV gives an example 
of the application of the theory. Finally, we summarize 
and conclude in Section V. 



II. WAVE EQUATIONS AND WAVEPACKET 
DYNAMICS: 



Consider light propagating in an inhomogeneous pho- 
tonic crystal characterized by the dielectric permittiv- 
ity tensor e (r) and magnetic permeability tensor /i (r) , 
where e and are spatially fast-varying periodic func- 
tions with slowly spatial modulation, e and ^ may also 
have smooth temporal modulation, which, together with 
dissipation of the medium, will be discussed in Sec. III.D 
for simplicity of the presentation. The Maxwell equations 
without external sources (i.e., currents) can be written 
as a Schrodinger-equation like form 



[H(-i9r,r)~S(c^,r)]$(r) = 0, 



where H ~ 



$ = 



idrX 
-idrX 
the electromagnetic wavef unction, 



(E,H)^ 
E{u) 



(1) 
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and the frequency lj corresponds 



( e{Lo) 

to idt in the time-dependent Maxwell equation. In mul- 
tiferroic materials with strong magnetoelectric coupling, 
S (cj) may also have non-diagonal terms (15j . The form 
of the wave equation ([T]) is very general. For instance, 
it can be taken as the Schrodinger equation for S = fto; 
and-H = -h^d^/2m + V {y). 

Without considering dissipation and frequency depen- 
dence of the medium, the wave equation ([T]) can be 
treated using the wavepacket dynamics. The central idea 
[3| is to construct a wavepacket 



1$) = dV(q,i)l^(q,rc)) 



with a mean wavevector 

qc = y d^(7|a(q,i)|^q 

and a preassigned wavepacket center position 
r, = [($|f|eo$) + (eo$|f|$)]/2p. 

Here ^(q,re) = {^,iif /^p - e^i-*'0(q,r, 
Bloch eigenstate of the local wave equation 

["Ho - So (wco, rc)]-0 (q, rc) = 



(2) 

(3) 

(4) 
is the 

(5) 



with the eigenenergy Wco — (q, Tc), — is the 
local operator with the required periodicity of the un- 
perturbed crystal, 8o = dEo/dcjco, 0(q, i"c) is the peri- 
odic part of the Bloch wave and satisfies {(f)\ Qol'P) = 
1, p = ((f>| 8o 1$) is the normalization factor. In a 
medium without frequency dependence, we have p = 
(E • D -I- H • B) /2, the energy density of photons. While 
in the Schrodinger equation, p is the wavedensity. The 
amplitude a (q, t) satisfy the normalization condition 

/A|a(q,i)|' = l. 

Using a time-dependent variational principle with the 
Lagrangian 



we obtain the semiclassical EOM for photons 
oqc 

,1.. 

rr'- c 



dVr 



(6) 

(7) 
(8) 



by following a similar derivation as that in Ref. 3 for 
electrons, where 



ujc = t^co + Im ( — 



{T-Lq — So) 



is the total energy, 



qq 



&0 



dqp/ 



c.c. 



is the Berry curvature in the momentum space, fiqr 
flrr have similar definitions. 



(10) 
and 



III. GEOMETRIC OPTICS APPROACH 

The geometric optics approximation we adopt does 
not involve the concept of wavepacket. Substituting the 
WKB trial wavefunction $ (r, t) = exp [iS (r, t)] u (r) into 
the time-dependent Maxwell equation, we obtain 



[H i-idr + p) - S (w, r)] u = 0, 



(11) 



where p = d^S is the wavevector, the frequency w = 
—dtS is the same as that defined in Eq. ([1]). p plays 
a similar role as qc in the wavepacket dynamics. In the 
geometric optics theory for an inhomogeneous medium, 
p's direction is perpendicular to the wave front that has 
a constant phase S. 

Because the spatial modulation of the medium is slowly 
varying, u can be taken as a periodic function of r over 
a few lattice sites. The spatial modulation of u can be 
accounted by taking u as a function of a slowly vary- 
ing variable R, where R describes variations on a length 
scale much larger than the lattice spacing. R plays the 
same role as Tc- For a fixed R, u is a periodic function 
of the fast varying variable r. With the slowly varying 
variable R and the fast varying variable r, Eq. (|lip can 
be rewritten as 



[H {~idf — idu. + p) — ^ {uj, r, R)] u = 0. 



(12) 



By choosing a suitable phase 5*, we can take p and uj as 
functions of R. 

Because u varies slowly with respect to R, we may 
solve Eq. perturbatively. Up to the first or- 

der correction 9rw, we can take the gradient expansion 
n i-idr - i^R + p) w -Ho i~idf + p)+An with 



and ^ {uj, r, R) = So (^o) + ©o'^i- Here we denote 



eo = 



dujo 



(13) 



(14) 



and expand the eigenenergy and wavefunction to the first 
order: uj = ujq +uji, u — uq + ui. With these expansions, 
the zero order of Eq. (fT2l) gives 



(15) 



[Ho {-idf + p) - So] uq = 0, 



which is similar as the local wave equation ([5]) in the 
wavepacket dynamics. For a fixed R, Eq. (llSp yields a 
Bloch wavefunction uq = y/pf, where we choose the nor- 
malization condition (/| Oo |/) = / d^r/* (r) 60/ (f ) = 
1. Here and later in this article, an inner product corre- 
sponds to an integration over a unit lattice cell. / plays 
the same role as (jj in the wavepacket dynamics. We 
emphasize that the wavefunction / is a periodic func- 
tion of the fast varying variable r that originates from 
the periodicity of the photonic crystals. / describes 
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the periodic wave property of hght in photonic crys- 
tals that does not exist in traditional geometric optical 
approach for an inhomogeneous medium. The periodic 
wavefunction can lead to Eq. (jlSp . together with suit- 
able boundary conditions, yields the energy dispersion 
relation luq = luq (p, R). In general / may be a function 
of r, R, p and wq, depending on the way to record the 
wavefunction. In this article, we will take / as a function 
of r, p, R only by replacing ujq in uq as ujq (p, R) using 
the dispersion relation. 

The first order of the gradient expansion of Eq. is 



(16) 



where £q — Hq — Sg, G ~ —AH + Qq^jJi- Note that the 
semiclassical geometric optics approach is different from 
the well-known Luttinger-Kohn treatment Tg"], which is 
a full quantum mechanical perturbative approach to per- 
turbed periodic systems. In Luttinger-Kohn treatment, 
the wavefunction is expanded around the band minimum. 
While in the geometric optics approach, the gradient ex- 
pansion is around the local Bloch wavefunction. 



A. Energy correction and Berry pliase 



Multiplying each side of Eq. (|16p with (uqI, we ob- 
tain the first order energy correction wi = (/| AH \f). A 
straightforward but tedious evaluation of uji yields 



T / 



Co 



(17) 



Co 



The first term wiq 
as that derived from the wavepacket dynamics 
second term ujn, = — Im (/| 8o 



^\ is the same 



In the 



^\ - R . ^ 4_ 

dt p dt OR 

p • ^ since we assume no explicit time-dependence of 
the medium. Therefore along a path C, a Berry phase 
7 = Ar • dR -f Ap • dp is accumulated. Here Ar and 
Ap are Berry connections in position and momentum 

spaces, with Ar = — Im (/| 9o -A-p has a similar 

definition. 

The procedure to obtain the energy correction (flT)) also 
yields the continuity equation 



dp 



aR-J=0, 



where 



9p 



E X H 



(18) 



(19) 



is the Poynting vector (the energy current density) of 
light. Note that p does not vary with time without con- 
sidering the dissipation of the medium, i.e., ^ =0. At 
the zero order, Jq — jOVg with Vq = duJo/d'p as the zero 
order of the local velocity. In the Schrodinger equation, 
J is the current density of particles. 



B. Poynting vector field 

Substituting u = uq+ui into ([T9|. we find the Poynting 
vector can be rewritten as 



i9p 



f^ppP - fipRVo ) -f (9r X pm. (20) 



Here u = ujo + t^ia, P = Vq • ^, ^r = ^ + ^ • 



m 



is the POAM defined in Eq. ([2T|) . The Berry curvatures 
57pp and f2pR have the same definitions as that in the 
wavepacket dynamics (Eq. ([TU])) with the replacement of 
the notations {fc, qc, Wco, </>} — {R, p,a;o,/}. The first 
three terms in Eq. ((20|) describe the energy flow due to 
the translational motion of photons with a velocity same 
as Eq. ([7]) for the wavepacket dynamics. 

It is clear from Eq. (I^Hl) that the energy current density 
J cannot be simply determined by the translational mo- 
tion of photons. J contains an additional term 9r x pm, 
which is absent in Eq. ([7]) for the wavepacket dynamics. 
Note that 



I Id] 
m = ( -— 

2\ap 



X (So - Ho) 



dl 
dp 



(21) 



is the POAM. For non-interacting electrons in crystals 
(no frequency dependence), m is just the orbital mag- 
netic moment of Bloch electrons [4] . In the wave-packet 
treatment of electron dynamics, m originates from the fi- 
nite spread of the wavepacket. The wavepacket generally 
rotates about its center position, giving rise to an orbital 
magnetic moment m = — e/2 {W\ (r — Tc) x v \ W), where 
I W) is the wave packet and v is the velocity operator [iTj . 
Therefore the term 9r x pm corresponds to a local cur- 
rent density field arising from the self-rotation of photons 
{i.e., POAM) in photonic crystals. 

Note that both p = d^S, and w — —dtS are gauge 
dependent quantities, i.e., they depend on the choice of 
S. In the static natural choice of the gauge is to 

ensure that energy is the same for different positions, i.e., 
dnuj = 0. Under this gauge, we find p = — which is 
the same as Eq. ([5]) in the wavepacket dynamics up to the 
first order of the gradient expansion. It is easy to check 
all terms in Eq. ([20)1 are gauge independent. Therefore 
the Poynting vector J does not depend the choice of S. 
The equations for J = pR and p determine the geometric 
ray self-consistently. With these two equations, we can 
obtain the position and momentum parameters R, p at 
any time. 

The gauge invariant momentum and energy defined 
in the wavepacket dynamics can be related with p, uj 

through qc = p- Ar, uJc = cjo+^i+Im (/| Oq Un- 
der this transformation, we can show that the equations 
for R and p transfer to Eqs. (O and ([5]) derived from the 
wavepacket dynamics (except the term 9r x pm). Note 
that the term J^rrR in Eq. (|5]) (with R replacing Tc) 
can be rewritten as R x K, where B 
effective magnetic field for photons. 



^ai3j^RR is an 
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C. Reciprocal electric field 

It is well known that the term — fjppp in the equa- 
tion of motion for R can be rewritten as — p x with 
n — iea^^Sl^^ as the reciprocal magnetic field in the mo- 
mentum space Q . Interestingly, a reciprocal electric field 
may also exist in a frequency -dep endent medium. In a 
frequency-dependent medium [l8|, if the wavefunction / 
is chosen in a way that has explicit luq dependence, the 
derivative ^ (x can be R, p) becomes ^ -I- ^^^^^ and 
many new terms depending on the Berry curvature in 
the frequency space may appear in the Poynting vector 
([20|) . For instance, the second term in Eq. ([20]) can be 
rewritten as 



pp X [Jl + vo X T] . 



where 



\dujo 



en 



} +C.C 



(22) 



(23) 



is the reciprocal electric field in the momentum-frequency 
space. It is interesting that the form of T obtained from 
the simple substitution here agrees with that obtained 
in Ref. |14l | for quasiparticle dynamics in an interacting 



Fermi liquid based on the complicated Keldysh formal- 
ism, where the self-energy of the quasiparticles has an 
energy dependence. 



and Poynting vector field due to the dissipation of the 
medium. For an isotropic medium (e and /i are scalar) 
with H = 0, we have 



= 2ti;oInie (wq) /Ree (wq) 



(25) 



We see Ren/2 corresponds to the attenuation con- 
stant or absorption coefficient of a plane wave propagat- 
ing in an absorptive medium [l9| . where n is the index 
of refraction of the medium. 

The dissipation of the medium yields an additional 
term to the Poynting vector field (l20t 



pF = plm 



{f\r-'Qo 



dl 
dp 



-2(/|r 



dl 
dp 



(26) 



which is also absent in Eq. ([7]) for the wavepacket dy- 
namics. Clearly this term is gauge invariant since t^^ = 
2(/|r|/). For an isotropic medium without frequency 
dependence, this term is zero. However, it can become 
nonzero in an anisotropic medium (e.g., Sxx 7^ £yy 7^ £zz)- 
pF is similar as the spin torque dipole density in spin 
transport theory (20| . where the existence of the spin 
torque reflects the fact that spin is not conserved micro- 
scopically in systems with spin-orbit coupling. The term 
pF originates from the dissipation of the medium, which 
induces the inconservation of the density p (i.e., p ^ Q). 
We thus name pF as the dissipation dipole density. 



D. Slow time modulation and dissipation of the 
medium 



In the presence of slow time modulation of e and p,, 
characterized by a slow variable T (similar as R), addi- 
tional terms J^Tq and — f^TR appear in the equations of 
motion ([7]) (also in (UHl) and ([5]) respectively. The term 
—VLt-r. — — {At is the Berry connection in the 
T space) can be taken as an effective electric field for 
photons with (— Ar, At) serving as the gauge potential 
of the effective magnetic and electric fields. Note that 
the effective magnetic field B and electric field — I^tr are 
different from the real magnetic and electric fields (E, H) 
which are contained in the wavefunctions of photons. 

In previous discussion, we have neglected the dissi- 
pation of the medium for simplicity of the presenta- 
tion. We note that the dissipation of the medium can 
be incorporated by adding an anti-Hermitian part iV = 
(S — S^) /2 to the operator S = S -t- iF, here F is as- 
sumed to be a small perturbation to the Hermitian part 
S = (S -|- S^) /2. 60 in previous equations is now re- 
placed with 80 — dSo/duj. The continuity equation ^TE\\ 
becomes 



drP + Vn ■ J = -p/t, 



(24) 



where = 2(/|F|/) describes the dissipation of the 
medium. The right hand side of the continuity equa- 
tion ([M)) represents the change of the energy density 



E. Degenerate bands 

In the above discussion, we consider only a single non- 
degenerate energy band. The generalization of the theory 
to degenerate bands (e.g., left and right circular polariza- 
tion states of photons) is straightforward 2l| . All above 
formulas are still valid by replacing the local wavefunc- 
tion / with Xafa, where is the Bloch wavefunction 
at band a and Xa is the superposition coefficient that can 
be determined through 



{U\H\fp)-i{fa.\eo 



dfi 
dt 



XH- 



(27) 



In a uniform medium without frequency dependence, Eq. 
((27| can be simplified as = — • ^^Xp on a heli- 
cal basis, which is the same as that derived in Q using 
wavepacket dynamics. Eqs. ()20|26|27p and the equation 
for p are the new set of the EOM for light propagation 
in photonic crystals. 



IV. TRANSVERSE ENERGY CURRENT FROM 
POAM 

Finally, we show that the additional energy current 
9r X /om in J can be nonzero and may significantly mod- 
ify the photon transport in photonic crystals using a con- 
crete example. Consider a 2D hexagonal photonic crys- 
tal with a uniform isotropic permeability tensor po^ij: 
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50 x/|G| 100 



Figure 1. (Color online) Plot of the trajectories of a photon 
in a 2D {xy) hexagonal photonic crystal obtained from the 
geometric optics approach. The length unit is \G\~^ = a/27r, 
where a is the lattice constant. A = 0.5, r;o = 0.008, drji/dx = 
4 X 10"'' with 771 {t = 0) = 0. The inital wavevector of the 
photon has p = (A/3,0) from the BZ corner K. Solid line: 
real trajectory of the photon. Dashed line: trajectory without 
dn. X jom. Dotted line: trajectory without both dn x pm and 
^prR. 



8(p2 + A2)2 dx ^ 
along the y direction. However, the nonzero m also yields 
another transverse energy current 



J'l = 9r X m 



Vic\G\ 



A2) drji 



8(p2 



dx 



(33) 



which is absent in previous literature. Note that these 
two currents have the same directions around p = and 
are comparable in magnitudes. In Fig. [l] we plot the tra- 
jectories of a photon with and without taking account of 
the contribution ([55)1 from the POAM. We see the clear 
distinction between these two trajectories, which indi- 
cates that the transverse energy current from the POAM 
cannot be neglected in the transport of photons in pho- 
tonic crystals. Finally we emphasize that the photon 
trajectories in Fig. 1 are obtained from our first order 
geometric optics treatment and it would be interesting to 
compare them with the fully correct photon trajectories 
obtained from ab initio calculation in the future. 



and an isotropic but spatially-varying permittivity ten- 
sor £ (r) eo, with [l^ 

(r) = £ (1 + <;Vg (r)) , (28) 
£xy (r) = ~£yx (r) = -i£ {vo + Vi^G (r)) • 

Here Vb (r) = 2 ^^^-^ cos (G„ • r), G„ are the three 
equal length reciprocal vectors in the xy plane, rotated 
120° relative to each other. The £xy (r) are the Faraday 
terms which explicitly break the time-reversal symme- 
try and lead to nonzero Berry phase contribution to the 
photon transport. 

We focus on the decoupled TE set {E^, Ey,Hz} of the 
electromagnetic fields. Around the corner K of the first 
BZ (see the inset of Fig. [T|), the photon dynamics are 
governed by an effective Hamiltonian [Toj 

H ^ uqI + c{pxax + PyCTy + Aa^) /2, (29) 

as shown in Ref. [l3|, where wq = c|G| (1 — <j/4) /VS, 
A = |G| — 3<^77o) /a/3, p is the deviation of the 2D 
Bloch vector from the BZ corner K, and c is the speed of 
light. From the Hamiltonian, we find 

n, = ±lA{p^ + A'y'^\ (30) 

m, = ±^A(p2 + A2)-\ (31) 

where ± correspond to two photonic bands. 

We assume the parameter rji (x) is slowly varying along 
the X direction. Consider a photon frequency at the lower 
band, a nonzero leads to a transverse energy current 

= — Oppp = — p X r2 



V. CONCLUSION 

In summary, we derive a new set of semiclassical EOM 
of particles in a chiral and dissipative periodic medium 
subject to slowly varying perturbations in space and 
time. The effects of several important properties of par- 
ticles (or waves) and media that are neglected in current 
literature {e.g., self rotation of particles, frequency (or en- 
ergy) dependence and dissipation of the media) have been 
included in the EOM. The EOM are derived based on a 
geometric optics approach for light propagation in peri- 
odic photonic crystals. We show that Berry curvatures 
in position, time, momentum, and frequency spaces, the 
dissipation dipole density of the medium, and the pho- 
ton angular momentum may affect the Poynting vector 
field of photons that traces the geometric optics rays. 
Reciprocal as well as effective electromagnetic fields for 
photons may arise from the Berry phase effects in pho- 
tonic crystals. We emphasize that although the geomet- 
ric optics theory and the EOM are presented for the light 
propagation in photonic crystals, they can also be ap- 
plied to many other systems, to name a few [ll|, electron 
dynamics in Bloch bands. X-ray propagation in natu- 
ral solid crystals that can be taken as photonic crystals 
for X-rays, light propagation in multiferroic materials, 
quasiparticle dynamics in a superfluid, and transport of 
acoustic waves. 
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